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Using 1 fb _1 of data from pp collisions at ^/s = 1.96 TeV at the Fermilab Tevatron collider 
collected by the DO detector, we search for decays of Kaluza- Klein excitations of the graviton in the 
Randall- Sundrum model of extra dimensions to e + e~ and 77. We set 95% confidence level upper 
limits on the production cross section times branching fraction which translate into lower limits 
on the mass of the lightest excitation between 300 and 900 GeV for values of the coupling k/M P i 
between 0.01 and 0.1. 
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PACS numbers: 13.85.Rm, 04.50.+h, 12.60.-i 

The large difference between the Planck scale, Mpi w 
10 16 TeV, and the weak scale presents a strong indication 
that the standard model is incomplete. In the presence 
of this hierarchy of scales it is not possible to stabilize the 
Higgs boson mass at the low values required by experi- 
mental data without an excessive amount of fine-tuning 
unless there is some, as yet unknown, physics at the TeV 
scale. 

Randall and Sundrum have suggested a model [l[ in 
which the fundamental scale of gravity is near the weak 
scale and gravity appears so feeble because it is expo- 
nentially suppressed by the existence of a fifth dimen- 
sion and a warped space-time metric. Standard model 
fields would be confined to one 3-brane (a 4-dimensional 
subspace of this 5-dimcnsional space) and gravity origi- 
nates at another 3-brane. Only gravitons propagate in 
the bulk between these two branes. The apparent weak- 
ness of gravity originates from the small overlap of the 
graviton wave function with the standard model fields in 
the fifth dimension. 

This model predicts a tower of Kaluza-Klein excita- 
tions as the 4-dimensional manifestation of the graviton 
propagating in 5-dimensional space. In the following we 
refer to these as RS (Randall-Sundrum) gravitons. The 
massless zero-mode couples with gravitational strength. 
The massive modes couple with similar strength as the 
weak interaction. Their properties are quantified by 
two parameters, the mass of the first massive excitation 
Mi and the dimcnsionlcss coupling constant to standard 
model fields, k/Mpi, where M pi = Mpi/\/8tt is the re- 
duced Planck scale. To address the hierarchy problem 
without the need for fine-tuning Mi should be in the 
TeV range and 0.01 < k/M Pt < 0.1 0. For these values 
the first massive RS graviton G is a narrow resonance 
with a width much smaller than the resolution of the DO 
detector. If kincmatically accessible, RS gravitons can 
be resonantly produced in high energy particle collisions. 
They decay into pairs of fermions or bosons. 

In this Letter we consider decays into e + e~ and 77 
pairs. We search for these as resonances in the e + e _ and 
77 invariant mass spectrum from 1 fb" 1 of data collected 
using the DO detector at the Fermilab Tevatron collider 
between October 2002 and February 2006. In the Teva- 
tron protons and antiprotons collide at y/s = 1.96 TeV. 
DO has previously published searches for RS gravitons [3j 
and excluded M x < 250_GeV for k/M Pt = 0.01 and 
M x < 785 GeV for k/M Pi = 0.1 at 95% confidence 
level with 260 pb" 1 of data. CDF has recently submitted 
for publication searches that exclude Mi < 889 GeV for 
k I'M pi = 0.1 based on 1.3 ftr 1 of data. 

The DO detector [j| @] consists of tracking detectors, 
calorimeters, and a muon spectrometer. The tracker em- 
ploys silicon microstrips close to the beam and concen- 
tric cylinders of scintillating fibers in a 2 T axial mag- 
netic field. The liquid- argon/uranium sampling calorime- 



ter has an electromagnetic section that is 20 radiation 
lengths deep, backed up by a hadronic section. The 
calorimeter is divided into a central section covering 
\r)\ < 1.1 and two endcap calorimeters extending cov- 
erage to |?7| < 4.2. The luminosity is monitored by two 
arrays of plastic scintillation counters located on the in- 
side faces of the endcap calorimeters. The pseudorapidity 
7] = — ln[tan(#/2)] and 9 is the polar angle with respect 
to the proton beam direction. The azimuthal angle is 
denoted by <fi and we measure object separation in the 
detector in terms of AR = \J (Ac/)) 2 + (Arj) 2 . We de- 
note the momentum component transverse to the beam 
direction with pr- Readout is controlled by a three- level 
trigger system. 

Since both electrons and photons result in electromag- 
netic showers with very similar signatures in our detector, 
we can define an inclusive selection that provides good 
efficiency for selecting e + e _ and 77 final states. In par- 
ticular we require clusters of energy depositions in the 
electromagnetic calorimeter that are consistent with the 
expected shower profile using a x 2 test and have less than 
3% of their energy leaking into the hadronic calorimeter 
section. We require that the cluster is well isolated with 
less than 7% of the cluster energy in an annular isola- 
tion cone with 0.2 < AR < 0.4 around the cluster cen- 
troid and less than 2 GeV for the sum of the pr of all 
tracks with 0.05 < AR < 0.4 with respect to the clus- 
ter centroid. To accept both electrons and photons we 
do not require a matched track. We start with a data 
set of 34 million events triggered on one or two electro- 
magnetic showers with px thresholds between 15 and 35 
GeV. We select events in which there are at least two 
such clusters with px > 25 GeV in the central calorimeter 
with 1 77 1 < 1.1. Including clusters in the end calorimeters 
would add little acceptance for decay products of massive 
objects. In the collider data we find 43639 events that 
satisfy these selection criteria with the invariant mass of 
the two clusters M ee / 77 > 60 GeV. 

Within the standard model, the Drell-Yan process and 
diphoton production give rise to e + e~ and 77 final states. 
The invariant mass spectrum for these is expected to fall 
towards higher masses except for the Z — > e + e~ reso- 
nance. We model these backgrounds using a Monte Carlo 
simulation with the PYTHIA [3] event generator using the 
CTEQ6L parton distribution functions § , followed by a 
GEANT-based @ detector simulation. Another source of 
events is the misidentification of one or two jets as elec- 
tron or photon candidates. The shape of the invariant 
mass spectrum of this source of events is estimated from 
data by selecting events with energy clusters in the elec- 
tromagnetic calorimeter that are not consistent with elec- 
tromagnetic showers and fail the % 2 test for the shower 
profile. The absence of the Z resonance in the back- 
ground spectrum in Fig. [T] confirms that this sample has 
no significant contamination from e + e~ final states. 
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We fit the shape of the invariant mass spectrum from 
the data near the Z resonance (60 < M ee / 77 < 140 GeV) 
with a superposition of the spectrum from Monte Carlo 
predictions for the standard model processes and the 
spectrum expected from misidentified clusters. In the 
fit, the spectra from e + e~ and 77 final states are nor- 
malized relative to each other by the leading order cross 
section from PYTHIA, the total number of events is fixed 
to the number of events observed in the data, and the 
fraction / of all events that have misidentified clusters is 
the only free parameter. We obtain best agreement with 
the data for / = 0.21 ± 0.01. The spectra are shown in 
Figure Q] Trigger thresholds affect the shapes near the 
low mass end of the fit window. We account for this by 
assigning a systematic uncertainty on the value of /. At 
masses above 100 GeV the trigger is fully efficient. 
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FIG. 1: Invariant mass spectrum from data (points). Super- 
imposed is the fitted total background shape from standard 
model processes including events with misidentified clusters 
(open histogram) and the fitted contribution from events with 
misidentified clusters alone (shaded histogram). 

We compare the invariant mass spectrum of our back- 
ground model with the fitted value of / to the data at 
higher masses. As shown in Figure [U we find agreement 
between background model and data in the high-mass 
range. There is a slight mismatch in the mass resolu- 
tion at the Z peak between our Monte Carlo simulation 
and the data. We verified that this does not affect the 
predictions of the background model at higher masses. 

From the fitted number of pp — > e + e~ + X events 
(most of them in the Z resonance), the acceptance and 
efficiency from the Monte Carlo simulation, and the cal- 
culated standard model cross section, we determine the 
integrated luminosity of the data sample. All Monte- 
Carlo derived efficiencies are multiplied by 0.96 so that 
the efficiency from the Z — > e + e~ Monte Carlo simula- 
tion agrees with the efficiencies measured in Z — > e + e~ 
data. The leading order cross section for the e + e~ final 
state with 60 < M ee < 130 GeV from pythia is 178 
pb. We multiply this by a next-to-leading order (NLO) 
iGfactor of 1.34 0. This gives 985 ± 35 pb" 1 . The un- 
certainty in this number is dominated by the uncertainty 
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FIG. 2: Invariant mass spectrum from data (points). Super- 
imposed is the fitted total background shape from standard 
model processes including events with misidentified clusters 
(open histogram) and the fitted contribution from events with 
misidentified clusters alone (shaded histogram). The grey 
shaded histogram shows the signal expected from gravitons 
with Mi = 300, 600, and 900 GeV and fc/M Pi = 0.1 on top 
of the total background. 



in the cross section from parton distribution functions. 
We do not include uncertainties on efficiencies and ac- 
ceptances because these cancel in the limit calculation. 
This value is in agreement with the number determined 
using the luminosity counters (1036 ± 63 pb -1 ) [Til ] . 

We determine the signal acceptance and efficiency us- 
ing a Monte Carlo simulation of RS gravitons with 200 < 
Mi < 1000 GeV using PYTHIA and GEANT. Systematic 
uncertainties in the signal efficiency originate from de- 
tector resolution (1-11%), parton distribution functions 
(0.2-5.5%), electron and photon identification efficiencies 
(1.4%), and the finite signal Monte Carlo sample size 
(0.5%). Contributions to the uncertainty in the back- 
ground prediction are from the finite size of Monte Carlo 
and data samples (2-24%), parton distribution functions 
(2-10%), the mass dependence of the NLO iGfactor (5%), 
and the uncertainty in the trigger thresholds (1%). In 
some cases the uncertainties vary with the invariant mass 
value. 

We compute upper limits for the production cross sec- 
tion of RS gravitons times branching fraction into e + e~ 
final states at 95% confidence level by comparing the ob- 
served and expected numbers of events in a sliding mass 
window. The width of the window was optimized for 
maximum sensitivity using the Monte Carlo simulation 
and varies from 20 GeV for M x = 200 GeV to 120 GeV 
for Mi = 950 GeV. We use a Bayesian approach to inte- 
grate over all important input parameters such as signal 
efficiency, background prediction, and integrated lumi- 
nosity, using a Gaussian prior with width equal to the 
estimated uncertainties in the parameters [12j. For the 
RS graviton production cross section we use a flat prior. 
To compute the limits, we use the integrated luminosity 
determined from the Z signal, which gives us a more pre- 
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FIG. 3: 95% confidence level upper limit on a(pp —* G + 
X) x B(G — > e + e~) from 1 fb _1 of data compared with the 
expected limit and the theoretical predictions for different 
couplings k/Mpi. 



cise normalization than the direct luminosity measure- 
ment. Figure [3] shows the limits as a function of invariant 
mass compared to predictions from the Randall-Sundrum 
model and Table Q] tabulates the results. Based on the 
observed and expected numbers of events we obtain lim- 
its on cr(pp — > G + X) x B(G — > e+e~/77k We divide 
by B{G -> e+e~/jj)/B(G -> e+e") = 3 13] to convert 
these to the quoted limits on a(pp G + X) x B(G — > 
e+e-). 

Using the cross section predictions from the Randall- 
Sundrum model with the same iVfactor as for the stan- 
dard model processes [HI , we set upper limits on the cou- 
pling k/Mpi as a function of Mi. This is shown in Fig- 
ure Hand tabulated in Table D For fc/M Pi = 0.01(0.1) 
we can exclude masses below 300(900) GeV at 95% con- 
fidence level. 

In summary, we have searched for RS gravitons as reso- 
nances in the e + e~ and 77 invariant mass spectrum from 
about 1 fb _1 of data from the Fermilab Tevatron collider. 
We find good agreement of the observed spectrum with 
standard model predictions and set lower limits on the 
mass of the first massive RS graviton at 95% confidence 



level of 300 GeV for k/M Pl = 0.01 and of 900 GeV for 
k/Mpi = 0.1. These are the tightest direct limits on RS 
gravitons to date. 
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